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Abstract. Geometric quantum discord, proposed by Dakic, Vedral, and Brukner 
■ [Phys. Rev. Lett. 105 (2010) 190502], is an important measure for bipartite 

correlations. In this paper, we generalize it to multipartite states, we call the 
(— | ' generalized version geometric global quantum discord (GGQD). We characterize 

Q_c GGQD in different ways, give a lower bound for GGQD, and provide some special 

states which allow analytical GGQD. 

, CT * PACS numbers: 03.65.Ud, 03.67.Mn, 03.65.Aa 

<N . 

>• . 1. Introduction 

o 

^ ' Quantum correlation is one of the most striking features in quantum theory. 

Entanglement is by far the most famous and best studied kind of quantum correlation, 
and leads to powerful applications [lj. Another kind of quantum correlation, called 
quantum discord, captures more correlations than entanglement in the sense that 
separable states may also possess nonzero quantum discord. Quantum discord has 
been attracted much attention in recent years, due to its theoretical interest to quantum 
theory, and also due to its potential applications [2 J . Up to now, the studies on quantum 
correlations, like entanglement and quantum discord, are mainly focused on the bipartite 
case. 

Quantifying the multipartite correlations is a fundamental and very intractable 
question. The direct idea is that we can properly generalize the quantifiers of bipartite 
correlations to the case of multipartite correlations [3j HI El EJ [7] . Recently, generalizing 
the quantum discord of bipartite states to multipartite states has been discussed in 
different ways |8j |9j LTOl EH]- As an important measure of bipartite correlations, the 
geometric quantum discord, proposed in [12J, has been extensively studied j2]. In this 
paper, we generalize the geometric quantum discord to multipartite states. 

This paper is organized as follows. In Sec. 2, we review the definition of geometric 
quantum discord for bipartite states. In Sec. 3, we give the definition of geometric global 
quantum discord (GGQD) for multipartite states, and give two equivalent expressions 
for GGQD. In Sec. 4, we provide a lower bound for GGQD by using the high order 
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singular value decomposition of tensors. In Sec. 5, we obtain the analytical expressions 
of GGQD for three classes of states. Sec. 6 is a brief summary. 

2. Geometric quantum discord of bipartite states 

The original quantum discord was defined for bipartite systems over all projective 
measurements performing only on one subsystem [131, [TJ] . That is, the quantum discord 
(with respect to A) of a bipartite state pab of the composite system AB (we suppose 
dimA = tia < oo, dim I? = ub < oo) was defined as 

D a (p A b) = S(p A ) - S(pab) + imn{S(U A (p A B)) ~ S(U A (p A ))}, (1) 

In Eq.(l), £(■) is Von Neumann entropy, pa = trsPAB, is a projective measurement 
performing on A, H-a(pab) is the abbreviation of (JIa®Ib)(pab) without any confusion, 
here Is is the identity operator of system B. Note that H-a^b^Pab)] = £?\b[IIa(pab)], 
that is, taking partial trace and performing local projective measurement can exchange 
the ordering. 

It can be proved that 

Da(pab) > 0, (2) 

riA 

Da(pab) = Pab = ^2pi\i)(i\® pf , (3) 

i=l 

where, ua = dimA, is any orthonormal basis of system A, {pf}i=i are density 

operators of system B, pi > 0, Yn=iPi — 1- 

The original definition of quantum discord in Eq.(l) is hard to calculate, even for 2- 
qubit case, by far we only know a small class of states which allow analytical expressions 

0- 

Dakic, Vedral, and Brukner proposed the geometric quantum discord, as [12] 

Da(Pab) = mm{tr[( P AB - v A b) 2 } : D A (a AB ) = 0}. (4) 

& AB 

Obviously, 

D g a (pab) = ^ Da(pab) = 0. (5) 

For many cases Dj[(pab) is more easy to calculate than Da(pab) since Dj[(pab) 
avoided the complicated entropy function. For instance, D^(pab) allows analytical 
expressions for all 2-qubit states [12], and also for all 2 x d (2 < d < oo) states [15]. 

3. Geometric global quantum discord 

In [IT], the authors generalized the original definition of quantum discord to multipartite 
states, called global quantum discord (GQD). Consider an iV-partite (N > 2) system, 
each subsystem A^ (1 < k < N) corresponds Hilbert space Hk with dimHk = n& (we 
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suppose rik < oo). The GQD of an iV-partite state p Al A 2 ...A N is defined as (here we use 
an equivalent expression for GQD [16J) 

N 

D(p Al A 2 ...A N ) = S (pA k ) ~ S(p AlA2 ... AN ) 
k=l 

N 

-moarnE S(U Ak (p Ak )) - S(U(p MA2 .,. AN ))}, (6) 

k=l 

where, II = H AiA2 ... An is a locally projective measurement on AiA 2 ...An. 
Similar to Eqs.(2, 3), we have Lemma 1 below. 
Lemma 1. 

D( PAi a 2 ...a n ) > 0, (7) 

D(p Al A 2 ...A N ) = p Al A 2 ...A N = Pii<a...ijvl*lX*l| ® l*2> <*2 1 ® •■■ ® |ijvX*iv|- (8) 

jii 2 ...ijv 

Where, { l*fc)}*l^=i is an y orthonormal basis of H k , k = 1,2,...N, r Pi Y i 2 ...i N > 0, 

Tliili 2 ...i N Pi\i 2 ...i N 1- 

Proof. Eq.(7) is proved in [11]. Eq.(8) can be proved as follows. Noting that 
^a 1 a 2 ...a n (pa 1 a 2 ...a n ) = U. Ai (U A2 ...(U An (p AiA2 ___ An ))...)), then by Eq.(3) and induction, 
Eq.(8) can be proved. 

With Lemma 1, in the same spirit of defining geometric quantum discord for 
bipartite states in Eq.(4), we now define the GGQD below. 

Definition 1. The GGQD of state p AlA . 2 ... AN is defined as 

D G (p Al A 2 ...A N ) = min {tr[p AlM ... AN - a AlM ... AN \ 2 : D{g Ai a 2 ...a n ) = 0}. (9) 

CTA 1 A 2 ...A N 

With this definition, it is obvious that 

d g ( PAiA2 ... An ) = ^ D{ PAiM ... An ) = 0. (10) 

In [17] . two equivalent expressions for Eq.(4) were given (Theorem 1 and Theorem 
2 in p2]), and they are very useful for simplifying the calculation of Eq.(4) and 
yielding lower bound of Eq.(4) [T71 [181 [19]. Inspired by this observation, we now derive 
the corresponding versions of these two equivalent expressions for GGQD. These are 
Theorem 1 and Theorem 2 below. 

Theorem 1. D g (pa 1 a 2 ...a n ) is defined as in Eq.(9), then 

D G (p Al A 2 ...A N ) = mm{tr[p AlA2 ___ AN - II(p AlA2 ... Ajv )] 2 } 

= tr[p AlA2 ... AN ] 2 - max{tr[n(p AlA2 ... j4iV )] 2 }, (11) 

where, II is any locally projective measurement performing on AiA 2 ...A^. 

Proof. In Eq.(9), for any o Ai a 2 ...a n satisfying D{(t AiA2 ... An ) = 0, o- Ai a 2 ...a n can be 
expressed in the form 

PA 1 A 2 ...A N = Puii--.wl l 'l)(*l| ® l*2> <*2 1 ® ••• ® Niv)(«iv|, (12) 

ili 2 ...iN 
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where, {|«k)}2 fc = i is any orthonormal basis of H k , k = 1,2, ...N. Pi 1 % 2 ...i N > 0, 
E lli2 ...i N Piii 2 ...i N = 1. We now expand p Al A 2 ...A N by the bases {|i fc )}" fc =1 = {b'fc)}"*=i, 
fe = 1, 2, ...iV. Then 

PAiA 2 ...Ajv = E Aiji,i2i2,..,WivKl)0'l| ® O2I ® ••• ® Niv)0'iv|, (13) 

«ljl,*2i2,---,j]vjjV 

fr|>A 1 A a „.A JV - O-Aaa-.A^] 2 = tr[{p AlA2 ... AN ) 2 ] + ~C feii 2 ...iiv) 2 



2 53 Phii,i2i2,---,iNiNPhi2 



. ,*2*2 , ■ ■ ■ ,«]V*]V-r «1«2 ■ ■ .IjV 

ilia— ijv 

^ r [(PAiAa— A;v) ] "I" 53 (Piiii,i 2 i2,— ,ijvijv — Piiia— ijv) 
iii 2 ...ijv 



- E (a 

lil,l3«3>— ,JiViiV/ - 

(14) 

iii 2 — ijy 

Hence, it is simple to see that when Pi 1 i 1 ,i 2 i 2 ,...,i N i N = Phi 2 ...i N for all ii, z 2 , ...,ijv, Eq.(14) 
achieves its minimum. 

Theorem 2. D g (pa 1 a 2 ...a n ) is defined as in Eq.(9), then 

D g (pa 1 a 2 ...a n ) = 53 



\2 

ia 2 ---ajv, 

aia2---ajv 



max 53 ( 53 A ai i 1 A a2 i 2 ...A aN i N C aia2 ... aN ) , (15) 

iii 2 ...i N a\a 2 ...a N 

where, C^^...^ and A* fc i fc are all real numbers, they are specified as follows. For any k, 
1 < k < N, let L(H k ) be the real Hilbert space consisting of all Hermite operators on 
iffc, with the inner product (X\X') = tr(XX') for X, X G L(H k ). For all fc, for given 

2 

orthonormal basis {X afc } a * ! of L(H k ) (there indeed exists such a basis, see [20]) and 
orthonormal basis {|«fc)}" fc =i of if^, C ili2 .. AN and are determined by 

PAiA 2 ...An 

E <?« 

X ai ®X a2 ® ... ®X ajv , (16) 

^ifc = (ifcl^ajifc). (17) 
Proof. According to Eq.(ll), and by Eqs.(16, 17), we have 

D G {p Al A 2 ...A N ) = tr[p Al A 2 ...A N } 2 - max{tr[n(p AlA2 ... yliV )] 2 } 
= E (c a ) 2 -max{tr[ ^ ^ C ai (ii\X ai \ii){i2\X a2 \i 2 ) 

a\a 2 ...a N iiia-.-iiv a i°2-c«)V 

. . . <ijv I-^ajv I * AT > I « 1> O'l I ® |*2)0"2| <S> ••• <8> |*Ar)(jjv|] 2 } 

) 1 (Caia2...ajv) HlclX ) ] ( ) ] ^aiii ^-a 2 i 2 •• •^ajviw^'aaaa ...ajv ) ■ (1^) 
aia 2 ---ajv ili2...i]v «l«2---"iv 



4. A lower bound of GGQD 

With the help of Theorem 2, we now provide a lower bound for GGQD. 

If we regard Pa 1 a 2 ...a n as a bipartite state in the partition {A k , A 1 A 2 ...A k _ 1 A k+1 ...A N }, 
then the original quantum discord and geometric quantum discord of Pa!A 2 ...a n with re- 
spect to the subsystem A k can be defined according to Eq.(l) and Eq.(4), we denote 
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them by D Ak {p AlM ___ AN ) and D% k (p AlA2 ... AN ). Comparing Eq.(3) and Eq.(8), it is easy 
to find that 

D g ( PAi a 2 ...a n ) =0^D% ( PAl A 2 ...A N ) = 0. (19) 
Consequently, comparing Eq.(4) and Eq.(9), we get 

D G (p AlA2 ...A N ) > D G Ak {p AlM ... AN ). (20) 

To proceed further, we need a mathematical fact, called high order singular value 
decomposition for tensors. We state it as Lemma 2. 

Lemma 2. [21 J High order singular value decomposition for tensors. For any 
tensor T = {T /3i/32 ^ 3n : f3 k G {1, 2, m k }, k = 1,2, N}, there exist unitary matrices 

= (Up klk ), such that 

Tf3 ll 3 2 ...f3 N = ^!7i^272'--^7iV A 7i72...7iV) ( 21 ) 

7i72---7iV 

V A* A = (22) 

i Sl72---7fc-i7fc7fc+i--7iV 7i72---7fc-l£fc7fc+l..-7iV i, 7 fc u 7fc£fc' K^) 

7172— 7fc-l7fe+l— 7iV 

4 fc) > > > > 0. (23) 
Combining Lemma 2, Eq.(20) and the lower bound of D A (p AlA2 ,„ AN ) in [17], we 

can readily obtain a lower bound of D g (pa 1 a 2 ...a n )- 

Theorem 3. D g (pa 1 a 2 ...a n ) is defined as in Eq.(9), then a lower bound of 

D G {p Al A 2 ...A N ) is 

tr[( PAl A 2 ...A N ) 2 } ~ min{ g s « : k = 1,2, N}, (24) 

7fc=l 

where are obtained by Lemma 2 in which let T = {C aia2 ,,, aN : £ {1, 2, = 
1,2, A/"}, C aia2 ... QiV are defined in Theorem 2. 

Proof. Since D g {j> AiA2 _ An ) and -D^.(paiA 2 ....4 ]V ) keep invariant under locally unitary 
transformation, hence the state Pa 1 a 2 ...a n in Eq.(16) and the state 

& AiA . 2 ... An = A 

X ai ®X a2 ®...®X aN , (25) 

have the same GGQD, and 

D A k (PA 1 A 2 ...A N ) = D Ak (A AlA2 _ AN ). (26) 
From Eq.(20), we have 

D g (A AiA2 ... An ) > D G k (A AlA2 ... AN ). (27) 
From the definition of D A (A Ai a 2 ...a n ), Lemma 2 and Theorem 1 in [T7J, we have 



D A k ( A AiA 2 ...A N ) = tr[(p AlA2 ... AN ) 2 } -max^( £ K^K^ 



Ha, 



) 2 

a 2 ...a N J 



= tr[(p AlA2 ... AN ) 2 ] - max £ A* < 

fife 

> tr[( PAlA2 ... AN r] - £ (28) 
We then attain Theorem 3. 



Geometric global quantum discord 6 
5. Examples 

We provide some special states which possess analytical GGQD. 
Example 1. For iV-qubit (N > 2) Werner-GHZ state 

p=(l-r)— + M)(i/,\, (29) 
the GGQD of p is 

D G (p) = p 2 /2. (30) 

In Eq.(29), I is 2 x 2 identity operator, p G [0, 1], is the iV-qubit GHZ state 

|V>) = (|00...0) + |11...1))/V2. (31) 

Proof. We prove Eq.(30) according to Eq.(ll). 
tr(p 2 ) can be directly calculated, that is 

Hp 2 ) = + P? + (2* - l)(^) 2 - (32) 

maxn{tr[n(p)] 2 } can be obtained by the similar calculations of Theorem 4 In [TB], the 
only difference is that the monotonicity of entropy function under majorization relation 
(Lemma 4 in [16]) will be replaced by the case of the function 

n 

f(pi,p 2 , -,Pn) = ~^Pi- (33) 
i=l 

That is, maxn{tr[II(p)] 2 } can be achieved by the eigenvalues 

r 1 -H , P 1 ~ P , V 1 z£_ 1 -/i 1 piA) 
2 N 2' 2 N 2' 2 N ' 2^ ' '"' 2 N 



Thus 

max{tr[II(p)] 2 } = 2(^ + ^ + (2 N - 2)(^) 2 . (35) 



Combine Eqs.(32, 35), we then proved Eq.(30). 
Example 2. For iV-qubit state 

p = 7^(i® N + cio-r + c^r + c 3 *ri m 



4 + 4 + 4 - max{c 2 , 4, 4} 
2 N 



the GGQD of p is 

D G ( P ) = ^ ' ^ ' ^ -p^l^l . (37) 

In Eq.(36), I is the 2x2 identity operator, {ci, c 2 , c 3 } are real numbers constrained by 
the condition that the eigenvalues of p must lie in [0, 1]. 

Proof. We prove Eq.(37) by using Eq.(ll). 

tr(p 2 ) can be directly found, that is 

Hp 2 ) = cl+ f N +cl - (38) 

maxn{£r[II(p)] 2 } can again be obtained similarly to Theorem 4 In [TB], the only 
difference is that the monotonicity of entropy function under majorization relation 
(Lemma 4 in [IB]) will be replaced by the case of the function Eq.(33). 
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Similar reduction shows max n {tr[Il(p)] 2 } can be achieved by each of them 

have multiplicity 2 Ar_1 , where c = max{|ci|, \c2\1 | C3 1 } . Therefore, 



max{tr[n(p)] 2 } = (39) 



1 + c 

"rr L "' L " vr/J J ~2 
Combine Eqs.(38, 39), we then get Eq.(37). 
Example 3. Anisotropic state 



P=0--*)-#T + s\<l>H<l>\, (40) 
the GGQD of p is 

D G {p) = s\l-\), (41) 

where, d = dimH 1 , H 1 = H 2 = ... = H N , I is the d x d identity operator, s G [0, 1], 
|0) = ^X)j=i \U---l), {\l)}f =1 is an fixed orthonormal basis of Hi. 

Proof. We prove Eq.(41) according to Eq.(9). For any locally projective 
measurement n, which corresponds iV orthonormal bases of Hi, we denote them by 
{|i fc )}? fc=1 ,A: = l,2,...,iV. Let {\l)}f =1 = {\m)} d m= i. Then, 

n(l0><0|) = - (ii\l)(m\ii)...(i N \l)(m\i N )\ii)(ii\® \i N )(i N \. (42) 

trWM - n(|0)(0|)] 2 } = 1 - 1 (J2(ii\l)(m\ii)...(i N \l)(m\i N )f 

> 1 - j2 E. E(^IOHn)-(^IOH^) = 1-5. («) 

and the minimum can be achieved by taking (ijfe|Z) = Si k) i, (m\ik) = S m> i h . We then 
proved Eq.(41). 

We make some remarks. For states in Eq.(29) and states in Eq.(36), the GQD can 
also be analytically obtained [TB], we then can compare the GQD and GGQD for these 
two classes of states. For states in Eq.(36) and states in Eq.(40), when N = 2, the 
GGQD in Eq.(37) and Eq.(41) recover the corresponding results in [22J. 

We also remark that, from Eq.(37), let the state Eq.(36) undergo a locally phase 
channel performing on any qubit, similar discussions as in [23] show that GGQD may 
also manifest the phenomena of sudden transition and freeze. 



6. Conclusion 



In summary, we generalized the geometric quantum discord of bipartite states to 
multipartite states, we call it geometric global quantum discord (GGQD). We gave 
different characterizations of GGQD which provided new insights for calculating GGQD. 
As demonstrations, we provided a lower bound for GGQD by using the high order 
singular value decomposition of tensors, and obtained the analytical expressions of 
GGQD for three classes of multipartite states. We also pointed out that GGQD can 
also manifest the phenomena of sudden transition and freeze. 
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Understanding and quantifying the multipartite correlations is a very challenging 
question, we hope that the GGQD proposed in this paper may provide a useful attempt 
for this issue. 
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